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ABSTRACT 


The presence of giant shells or loops in giant HII regions are clear witness of 
the mechanical energy input from massive stars. Here we evaluate the impact 
that winds may have on the structure of giant nebulae and on their supersonic 
velocity dispersion. We follow the suggestion from Chu & Kennicutt (1994) to 
see if a combination of a large number of unresolved wind-driven shells caused 
by massive stars could produce the integrated broad Gaussian profiles typical of 
giant HII regions. 

The results, accounting for a wide range of energies, densities and velocity 
or age of the expanding shells, show that supersonic Gaussian profiles may arise 
only from a collection of unresolved wind-driven shells if the shells present a 
peculiar velocity distribution which implies a strongly peaked age distribution 
leading to an awkward star formation history. On the other hand, a uniform 
distribution of ages originates profiles with a fiat-topped core defined by the 
terminal shell velocity and a steep decay as up to the largest detectable 
shell speed. Thus, supersonic profiles can arise only if the final speed of the 
unresolved shells is supersonic. This implies an equally supersonic random speed 
of motions in the ionized gas disrupting the shells before they slow down to 
subsonic speeds. It also implies a mechanism, independent of the shells caused 
by massive stars, responsible for the supersonic stirring of the background 
medium. These facts together with the conditions for shells to remain unresolved 
to present day devices (energies, final speeds and ages), indicate that the winds 
may be produced by low-mass stars. In this latter case, if the sources move 
supersonically in the gravitational potential of the whole system they could stir 
the gas, with their cometary bow shocks, to a Ugas ^ (^stars causing a supersonic 
local random speed of motions within the system. 
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1. Introduction 

One of the most debated issues in the field of giant HII regions is without doubt 
the nature of their supersonic gas motions. It is well understood that an energy source 
is necessary to replenish, at all times, the dissipation caused by radiative cooling behind 
the multitude of shock waves that should result from the supersonic motions. In this 
respect the integrated mechanical energy deposited by the large number (hundreds or 
thousands) of massive stars, in most cases inferred from photoionization requirements, has 
been presumed to be the main mechanism driving supersonic turbulence (Dyson 1979; Rosa 
& Solf 1984; Hippelein & Fried 1984; Chu & Kennicutt 1994, hereafter CK). Given the 
apparent association between molecular gas and the giant nebulae, champagne flows have 
also been thought to be responsible of the supersonic random motions (Skillman & Balick 
1984; Hunter & Gallagher 1985). 

On the other hand, the supersonic velocity dispersion (a) in giant HII regions has been 
shown to correlate with the size (a ^ size°-^) and luminosity (a ^ of the giant nebulae 

(Terlevich & Melnick 1981), and the scatter in the data is strongly reduced if one considers 
instead of the core or the halo dimensions, or a combination of the two, the ’’kinematic 
radius” (Munoz-Tunon 1994), i.e. the dimension of the area where high-resolution studies 
warrant a consistently uniform supersonic Gaussian profile independent of location. These 
facts have no apparent explanation within the framework of the energetics caused by 
massive stars. Terlevich & Melnick (1981) noted that the same correlations apply to 
virialized systems such as globular clusters and the nuclei of elliptical galaxies, and from this 
they suggested that gravitational energy must be the main source that drives supersonic gas 
motions, accounting in this way for the observed correlations. In their original suggestion 
Terlevich & Melnick proposed as a possibility a large number of ionized condensations 
moving supersonically in the virialized (gas + stars) system, and called it the ’’gravitational 



model”. The model however requires a mechanism to replenish the ionized fragments, given 
that champagne effects would disperse condensations in a time much shorter than the HII 
region life-time. Their observational resnlt, however, has been thoronghly confirmed and 
also extended to HII galaxies (see Melnick 1992 and references therein). 

Turbulence has also been suggested as a mechanism to explain both the supersonic 
width and the observed relationship between velocity dispersion and giant HII region size 
(Roy et al. 1986). There are however, several important issues regarding this suggestion. 
For instance, it is not clear from the velocity structure that the observed random motions 
are evidence of tnrbnlence, and, as observed by Castaneda (1988), one is not even snre 
that this ’’tnrbnlence” is the same as that described by the conventional theories of finid 
mechanics (see also Melnick et al. 1987). 

A more recent explanation of the giant HII region empirical correlations was derived 
as a direct consequence of the formation of massive stellar clusters ( Tenorio-Tagle , 
Munoz-Tunon & Cox 1993). There the supersonic velocity dispersion arises from the 
virialized motion of an ensemble of recently formed low-mass stars which nndergoing 
low-energy winds, prodnces a snfficient nnmber of cometary bow shocks to stir the 
surrounding matter to a agas ^ <7stars and thus locally enhance the random speed of motions 
to a snpersonic valne. The virialized passage of the cometary shocks is responsible for the 
stabilization of the collapse and also nltimately for the observed a vs size and Inminosity 
correlations found by Terlevich & Melnick (1981). In the case of massive clusters (M > 10® 
Mq), the motion of the virialized stars and of the random speed of motions the parent clond 
remain snpersonic (> 10 kms^^) even after massive star formation and photoionization. In 
the cometary stirring model, during this phase the parent cloud, formerly only subjected to 
the cometary passage of the many low-mass wind-driven sources becomes disrupted by the 
mechanical energy of the massive members of the clnster. And thns, as the volnme affected 
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by the mechanical energy from massive stars grows with time, matter is removed from the 
location where the cometary sources could stir it, finally losing the information about the 
true origin of the supersonic velocity dispersion. 

Recently, new data sets on 30 Doradus (CK) and NGC 604 (Sabalisck et al. 1995) 
have refreshed the interest on the topic. The NGC 604 2D Fabry-Perot data set spatial 
sampling is superior to any data set ever taken of 30 Doradus, including that of Smith & 
Weedman (1970) or that of CK, and it surely is useful to clarify the nature of the giant HII 
region supersonic gas motions. In the case of 30 Doradus the observations of CK revealed 
a number of major kinematic subsystems (shells, rings, filaments, networks of shells, etc.) 
and a line splitting of 20-150 kms^^. From this, the authors concluded that the most 
important physical mechanism which produces the global velocity dispersion is that which 
accelerates the shells, namely stellar winds from OB associations. The result from NGC 
604 seems totally different. There, all lines clearly arising from shells and or filaments, 
apart from being strongly asymmetric or even split lines, are all very low intensity lines. On 
the other hand, lines arising from the two brightest emitting knots in NGC 604 are broad 
well behaved Gaussians implying a supersonic velocity dispersion intrinsic to the brightest 
knots, which do not seem at all dominated by the action of the many strong stellar winds 
operating in the region. 

Here we evaluate the impact that winds may have on the structure of giant HII regions 
and on the supersonic velocity dispersion. Section 2 analyses the outcome of recent high 
resolution studies on the nearest giant HII regions. We follow the suggestion of CK looking 
to see whether the integrated profile caused by hundreds of unresolved wind-driven shells 
could combine to produce broad Gaussian profiles. Section 3 derives the emission profile 
produced by single and multiple shells, accounting for a wide range of powering energies, 
background densities and velocity or age range of the expanding structures. Section 4 lists 
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our conclusions. 

2. Supersonic turbulence in the nearest giant HII regions 

The following discussion is based on recent studies of the two closest giant HII regions, 
30 Doradns at a distance of 50 kpc in the LMC, and NGC 604 in M33 at a distance of 
720 kpc. The two objects have been the subject of recent high-resolution spectroscopy and 
the corresponding outputs are used here to draw firm conclusions regarding the supersonic 
natnre of the velocity dispersion in nearby giant HII regions. 

In the case of 30 Doradus, the inner 9’^, corresponding to the central 135 pc, were 
covered by more than 37 echelle spectra (CK). From these, many of the wind-driven shells 
and networks of shells prodnced by the mechanical inpnt from massive stars were clearly 
identified through emission-line patterns. This led the authors to conclude that “at high 
spatial and spectral resolution the nebular velocity field breaks up into hundreds of discrete 
snpersonic strnctnres with diameters ranging from less than 1 to 100 pc, expansion speeds 
of 20-300 kms^^, and energies of lO'^^^^^erg.” 

The main point in their stndy is however, that the snm of the 37 spectra from the 
core region (some of which arise from well identified shells or sections of shells) resnlts 
in a well behaved Gaussian profile (their figure 11). This led to the conclusion that “the 
unusually high supersonic velocity of 30 Doradus when measured on large scales is due to a 
combination of hnndreds of discrete motions combined with a smooth nnderlying tnrbnlent 
velocity field”, and by extrapolation to conclude that “the use of global velocity profiles of 
giant HII regions to draw inferences about the physical origin of the gas motions is a futile 
attempt.” 


There are, however, several points that seem to have been overlooked in the analysis 
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of CK, and that we regard as fundamental to this discussion. The background or ’’smooth 
turbulent velocity field” seems totally free or uncontaminated by the presence of any of the 
resolved shells, and thus one should notice that : 1) the so called smooth turbulent velocity 
field found as a background component in the halo, and with the same properties (line 
width) in the core of the region, also presents a Gaussian profile; 2) the derived a value 
of this component is also supersonic {^background ^ 18 kms^^), and 3) such a component 
accounts for at least half (some 5.6 10^^ erg) of the kinetic energy present in 30 Doradus. 
The smooth Gaussian profile of the background component is clearly seen in their figures 
9a and c, which correspond to typical spectra from the halo (far away from R136) and the 
central core, respectively. The two lines differ by a factor of about 400 in intensity but 
otherwise they are identical with a supersonic a value of 15-20 kms“^. These lines are 
the ones CK account for with a collection of unresolved wind-driven shells. Note on the 
other hand that, if supersonic gas motions are due to gravity, the measured a value of the 
background component implies a total gravitational mass Mxotai = 220{R/pc) {a/kms^^) 
(see Terlevich & Melnick 1981) in the range 3-6 xlO® Mq , if one uses R = 67.5 pc (the 
radius of the central 9’), and a a in the range 15-20 kms^\ which is consistent with the 
total detected mass of the system (Melnick 1992). 

From the 2D Fabry-Perot study on NGC 604 (which led to 65000 spectra in , 
and 65000 in [OIII] over the central 230 pc of the nebula). Sabalisck et al. (1995) derived 
the particular sectors of the nebula that cause the global velocity dispersion value of this 
source. It was shown that the brightest knots (knots 1 and 2) are the ones that show clearly 
well behaved Gaussian profiles with a a value (^ 17 kms^^) similar to the global velocity 
dispersion associated with NGC 604 ((T 604 ) by earlier single aperture studies. The two 
emitting knots (^ 7 pc in radius) lie within the central zone of NGC 604 separated by 50 pc 
from each other. Many other values of a were found in the region, both subsonic and even 
larger than (T 604 . The latter ones result from a Gaussian fit to strongly asymmetric or split 
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lines, most likely caused by the photoionized matter swept up by the winds from the several 
W-R and massive stars identified in the region. However, in all cases the intensity of those 
lines is much smaller than that presented by the truly Gaussian profiles that emanate from 
either knot 1 or knot 2. Clearly, in NGC 604 the sampling of the brightest zones leads to 
profiles similar to the integrated global ones. This set of observations implies the existence 
of a motor intrinsic to knots 1 and 2 capable of causing the same supersonic (T 604 value in 
both of them, despite their physical separation. For the case of 30 Doradus the background 
gas shows the same asoBor both in the core and the halo implying an extended energetic 
mechanism in continuous operation. CK inferred by extrapolation of their integrated line 
profile that the supersonic Gaussian lines emanating from the background flow should 
also be the result of a large number (hundreds) of spatially unresolved wind-driven shells 
from massive stars. Here we follow their suggestion, looking to see if a large number of 
unresolved shells does indeed lead to the supersonic Gaussian lines typical of the background 
component of giant HII regions, and whether this has any further implications on the 
history of star formation. Our approach differs from the original one of Dyson (1979), in 
that here different powering energies, background densities, and age distributions of the 
wind-driven shells are considered. 

3. The signature of unresolved multiple shells 

From the equations of motion and evolution of wind-driven shells one can derive the 
profiles produced by multiple shell systems. Here it is further assumed that the predicted 
line profiles result from a population of wind-driven shells acting concurrently on an HII 
region while remaining unresolved to the high spatial resolution of present-day devices. The 
major goal of this exercise is to test whether such a system is capable of producing the 
smooth, supersonic Gaussian lines observed in giant HII regions. 
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3.1. Individual shells 

The mechanical energy deposition from massive stars leads to a double shock pattern. 
The inner shock decelerates the wind while thermalizing its kinetic energy and the outer one 
sweeps up the surrounding background gas (see Weaver et al. 1977). The two gases are well 
separated by a contact discontinuity, and are kept in pressure equilibrium throughout the 
evolution. Thus, at the onset of strong radiative cooling, the matter swept up by the outer 
shock condenses into a geometrically thin outer shell. In the presence of a large ionizing 
photon flux, this will be kept fully ionized at a temperature ^ Thii, while contributing to 
an emission line spectrum typical of HII regions. Wind-driven shells around massive stars 
(or groups of stars) thus develop isothermal shocks, and consequently the density in the 
shells is times the ISM density, where A4 is the outer-shock Mach number. Under these 
conditions, the equations describing the evolution of the velocity (vg), outer radius {Rg), 
thickness (ARg), mass (Mg) and Ha luminosity (Lg) of a single isolated shell are {e.g. , 
Dyson 1979): 


Vg{t) 

Rg{t) 

ARg{t) 

Mg{t) 

Lg{t) 


100 (e/n)i/^t-2/^ km s^^ 
1 .6(e/n)i/®t^/® pc 
Rg{t) 




1.6 X 10^"‘(cij///10)^^nte erg s^^ 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where e is the mechanical luminosity of the wind in units of 10^® ergs“\ n is the ISM 
density around the star (in cm^^), t is the age of the shell in units of 10"^ yr and chu is the 
sound velocity (in kms“^). Eq. (4) assumes a mass per particle // = 0.6mp corresponding 
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to a fully ionized gas with solar abundance. The numerical coefficient in eq.(5) assumes a 
temperature of 10"^ K, and hence a sound speed chii ^ 10 kms“^and an Ha recombination 
coefficient of 8.7 x 10^^"^ cm^^ s^^ (Osterbrock 1989). The colnmn density throngh the shell 
is given by 


S, = n,AR, = 1.7 X - {M‘^ - 1)^/^] cm-^ (6) 

For n = 100 cm“^ and e = 1 the shell column density increases from 1.8 x 10^® to 
1.7 X 10^^ cm“^ as the shell slows down from Vs = 100 to Vs = 10 kms“^. Assnming that 
only abont 10^"^ — 10^^ of the H atoms are in the level n = 2 this yields snch low valnes of 
the Her optical depth as clearly warrant the optically thin approximation. 

The line profile prodneed by a constant velocity (wg), optically thin spherical shell is 
flat-topped between —Vs and Vs (e.g. Beals 1931, see Fignre 1). Thermal motions of the gas 
smooth the profile on scales of ~ chii, but otherwise the emission should be constant across 
the line profile, with an intensity given by Fs{v) = Ls/2vs (see Fignre 1). 

A number of effects could in principle alter this profile shape, (i) At large optical 
depths radiative transfer effects come into play and alter the shape of the line profile. 
However, the Ha optical depth is small, and so the optically thin approximation is a valid 
assnmption. (ii) If the radial velocity is not constant across the shell then a different profile 
would result. Such velocity differences are not expected in wind-driven shells given their 
small geometrical thickness, and thns small velocity gradients wonld not be snfficient to 
prodnee significant profile changes, (in) Non-spherical, partly broken or one-sided shells 
are often observed in giant HII regions. Sections of a spherical shell (“caps”) still produce a 
fiat top profile, but the velocity range would be smaller than —Vg Vg, and dependent on 
the orientation of the shell section. In a system containing many shells, however, one wonld 
expect a shell section oriented in one direction to be compensated by another shell section 
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oriented in the opposite direction. The final effect would be that the profile resulting from 
a number of incomplete shells is equivalent to the profile of a single complete shell. This 
equivalence should apply unless the overall symmetry of the system is not spherical, (iv) 
Large shells in nearby giant HII regions are often spatially resolved, and their line profiles 
show line-splitting unless the whole shell is included in the spectrograph aperture. 

3.2. Multiple shells 

The line profile resulting from the superposition of many shell profiles can be computed 
using the formalism developed by Cid Fernandes & Terlevich (1994) and Cid Fernandes 
(1995). In a system containing many wind-driven shells (Ng), the total line profile will 
be the sum of many fiat-top profiles, each with a different width and intensity, as these 
corresponds to different ages, mechanical luminosities and ISM densities: 


Ns 

Ftot{v) = Y, Fs{v; ti, Cj, Hi) (7) 

i=l 

For a large number of shells, this summation can be approximated by an integral over the 
distribution of t, e and n: 


tot 


/ Fs{v; t, e, n)dNg{t, e, n) = Ng /// Fg{v,t,e,n)p{t,e,n)dtdedn (8) 


where Ng is the total number of shells and pit, e, n) denotes the probability density of a 
given set of t, e and n. A large variety of profile shapes can be obtained by arbitrarily 
choosing the probability distribution function, and each of these (see below) results into a 
very precise shape of the total line profile. 
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3.2.1. Time-averaged line profiles 

The probability p{t,e,n) can be written as p{t\e,n)p{e,n); i.e. , the prodnct of the 
probability that a shell with a given e and n is present at a time t and the distribution 
function of e and n among all shells. Eq. (8) then becomes 


Ftot{v) = Ns J j Ftot{v;e,n)p{e,n)dedn (9) 

where Ftot{v, e,n), is the average prohle of shells with a given e and n but different ages, 
given by 


Ftot{v;e,n)= Fs{v;t, e,n)p{t\e,n)dt (10) 

Jti 

where t* and tf are the initial and hnal ages of the shell respectively and F is the age 
of a shell with velocity Vg = v (older, slower shells do not contribute to the prohle at 
v). A shell can only emit in Hct if it is subjected to photoionizing radiation. Thus, 
a necessary requirement is that the shock causing the shell slows down to velocities 
~ 100-250 kms“^or less to ensure rapid recombination. Faster shocks (younger shells) lead 
to higher temperatures in the swept up matter, and given the interstellar cooling law enter 
a quasi-adiabatic regime that inhibits recombination. This sets the initial age to t* ~ t(vs = 
100 kms“^) = 10"‘(e/n)^/^ yr. The hnal age tf should correspond to the time at which the 
shell velocity approaches the random speed of motions of the interstellar medium, losing its 
identity upon disruption (Spitzer 1968). 

Regardless of the exact values of and tf, replacing Fs by Ls/2vs in equation 10 and 
using eqs. 1 and 5 yields the following solution for Ftot{v, e, n): 
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Ftot{v;e,n)= ^ / p(t|e, n)dt = n)dt (11) 

hi ^Vs[t) Jti 

where A = 8.0 x 10^^(cij///10)“^ ergs^Y^ms"^. Therefore, the profile at velocity v is 
proportional to the mean of the shells faster than v. 

The exact distribution of ages is unpredictable without the detailed history of star 
formation in the system. Information on the time at which strong stellar winds begin to 
act for stars of different masses, and the span of time that they remain in full operation 
would also be required to evaluate p(t|e, n). However, unless star formation is synchronized 
to less than the life-time of a shell (r^ ~ 10® yr for massive stars), it seems a reasonable 
approximation to assume all ages to be equally likely, with old shells from the first stages 
of star formation and young shells from newly born stars occurring in similar proportions. 
The probability of observing a shell with age between t and t + dt is therefore p{t)dt = dt/r, 
where t = Tgf + Tg is the duration of the star-formation episode plus the lifetime of the 
shell. Both and can in principle depend on e and n. This distribution of ages yields 
the following solution to eq. (11) 


^ / X t m12/5 12/5i 

Ftot{v;e,n) = - — - tf)] - t] } 


( 12 ) 


which, making use of the relation between Vg and t (eq. 1), is equivalent to 


Ftot{v,e,n) = 


5Ae z' Vi \ ® 


12t VlOO 


maxr; 


1^/)' 


- 1 


(13) 


Therefore, the total line profile expected from systems containing many wind driven 
shells with an uniform distribution of ages has a fiat-top core between —Vf and Vf, and 
decays as up to the maximum detectable shell velocity, Vi (see Figure 2). This very 
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steep profile is a result of the rapid decay of Vg with time and the even more rapid increase 
in line luminosity as t increases, which combined give a much larger weight to older, slower 
shells. Note also that since all shells contribnte to the line core and old shells contribnte 
only to the line core, the total profile necessarily peaks at v = Vf. 

The hypothesis that all ages between t* and tf occur with equal probability in multiple 
systems plays an important role in shaping Ftot{v; e,n). Thongh this seems to be the 
most natural assumption, one might wonder what line profiles would result from different 
age distributions. In the extreme case where all the stars form simultaneously and all 
their winds act at the same time, all shells wonld be observed with the same age t, and 
Ftot{v; e,n) wonld be identical to Fs{v;t, e,n), i.e. a fiat-top profile between —Vs{t) and 
Vs{t). A power-law distribution of ages with p{t) oc would produce a profile similar to 
eq. (13), but with wings decaying as instead of 

One can also invert the problem and ask which distribution of ages/velocities would 
result from a chosen line profile. This is most easily done re-writing equation 11 as an 
integral over shell velocities instead of ages. Using eq. 1 one obtains 

Ftotiv] e, n) = [ v^'’l‘^p{vg)dvg 

JTasxivf ,\v\) 

which can be inverted to 

(15) 

avg 

This solution applies only for |u| > u/, i.e. outside the line core. The profile between 
—Vf and Vf cannot be changed, being necessarily fiat-topped irrespective of p(us). The 
distribntion of shell velocities reqnired to prodnce a Ganssian profile with standard 
deviation a is shown in Figure. 3, along with the corresponding distributions of ages and 


(14) 
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the resulting line profile. (Note that the profile is Gaussian only outside the —Vf to +Vf 
core, where it remains fiat-topped). In this case p{vs) oc exp{—v‘^/2a‘^), which peaks 
at Vs = (9/2)^/^(t = 2 .12(7. Thns, in order to obtain a Ganssian profile a very peaked 
distribntion of shell ages is reqnired. Snch an age distribntion implies a very contrived and 
hard to explain star-formation history. It could also imply a rather unlikely event given the 
short period of operation (< 10^ yr) compared to the life-time of giant HII regions (< 10^ 

yr)- 

In what follows we shall assume a uniform distribution of shell ages throughout the 
life-time of the nebulae, a possibility that seems a more natural distribution than the ones 
shown in Fignre 3c. 


3.2.2. Variable e and n 

Strictly speaking, eq. (13) applies only to a population of shells with identical driving 
energies and evolving in a nniform-density medinm. It is well known however, that 
wind-driven shells in giant HII regions have a variety of values of e and clearly evolve into a 
highly non-uniform density medium. These could in principle modify the line profile shape. 
Note however that the total profile after convolving eq. (13) with p(e, n) wonld change its 
shape only if the initial and final velocities Vi and Vf have a dependance on e or n. In 
fact, eq. (13) remains valid after e/r is replaced by the mean value of this ratio among the 
varions shells. If p(e, n) is to affect the shape of the total line profile it has to do so by 
affecting Vi or Vf, and thns one might as well replace p{e,n) by p{vi,Vf). Given that the 
profile Ftot{v;e,n) is so strongly concentrated around the terminal velocity, variations in Vi 
from shell to shell would not have a strong impact on the total profile, and thus one can 
fnrther simplify the problem by considering only variations in wy. With these considerations 
the total line profile can be written as 
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where and define the range of possible Vf valnes. To compnte the total profile one 
needs to specify how e and r depend on n/, besides the probability distribntion of n/ itself. 
Since these functions cannot be derived from first principles, a parameterized approach is 
required. Eq. 16 can be evaluated by defining a weight function w{vf) = e{vf)p{vf)/T{vf) 
and parameterizing w as a power law: w{vf) oc vj. The resulting line profile is fiat topped 
from —to and follows a v~^ law for |n| > vY^■ For velocities between and 
the profile is a combination of two power laws, one with a —6 slope and the other with a 
slope of 7 — 5. Figure 4 illustrates this profile for different combinations of 7, Uy”"’ and v^^- 

There are no strong conclusions to be derived from this analysis without a physical 
justification for the w{vf ) expression. Nevertheless, Figure 4 already shows that unless the 
weight function is very strongly biased towards large terminal velocities, the total profile 
shape remains far from Gaussian. Roundish profiles begin to arise for 7 > 7. Note, however, 
that even if one could physically justify such possible values of 7, a large dynamic range of 
Vf^s is also required to produce a total Gaussian-looking intrinsic profile. 


3.3. Unresolved shells 


For a given value of the ratio e/n, the constraint placed by spatial undetectability of 
wind-driven shells {Rg < Vcrit) in giant HII regions can also be written as a limiting velocity 
criterion: 
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= 31 km/s 

Clearly, low values of tjn {i.e. either low energy winds and/or high densities) favor the 
constraint as well as young, and thus small and/or fast moving, shells. The constraint is 
illustrated in Figure 5, where the minimum velocity (Vcrit) such that shells are not spatially 
resolved are indicated for a range of e/n. Shells below the lines drawn for different 
valnes wonld have exceeded Vcrit and thns be resolved. In the hgnre, the possible range of 
e/n values to account for the 30 Doradus “slow shells” resolved by CK is clearly indicated 
for an assnmed = lOpc. The range of valnes reqnired to acconnt for faster shells lies to 
the right and ont of the plot, thns implying the joint action of many massive stars powering 
the detected fast expanding structures. Shells above the lines are unresolved. Note however 
that the rapid wind-driven shell evolution implies that unavoidably, most of the unresolved 
shells will soon have velocities close to Vcrit- If one was to exclnde, as demanded by the 
constraint, resolved shells from the computation of the total line prohle arising from the 
background in giant HII regions, one could treat Vcrit as the terminal velocity vj. This 
together with the valnes of velocity dispersion (15-40 kms^^) and backgronnd density 
(1 < n < 10^ cm^^) fonnd in giant HII regions implies rather low valnes of e, in fact 
indicative of the power expected from low mass stars {E^ind ^ 10^^“^®ergs“^) as shown by 
the shaded area in Fignre 5. The above conclnsion becomes more stringent the smaller the 
valne of rcru- 


3.4. Synthetic observed profiles: thermal and instrumental effects 

Here we consider the fact that observed line profiles are a convolntion of intrinsic 
profiles with both the thermal motions of the photoionized gas and the instrnmental profile. 
Both thermal and instrumental effects can be represented by Gaussians with a combined 
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width given by (Tq = {(7th + Typical values of ath and atnst are 10 kms^\ which 

yield gg = 14 kms“^. In principle, the intrinsic profile could be recovered by applying 
a deconvolution algorithm to the observed profile, but that is seldom done. Instead, 
the usual practice is to present the observed profiles and quote their corrected width 
(7cor = under the assumption that all profiles considered behave as Gaussian 

distributions. Having this in mind, one can predict the shapes of observed profiles and their 
corrected width for different combinations of the intrinsic profile width and the expected 
values of gg- 

Figure 6 shows a mosaic of profiles for different values of u/ (which sets the width of 
the theoretical profile) and gg- The thin lines represent the theoretical intrinsic emission 
from a large number of unresolved shells, as calculated from equation 13. The dotted 
lines are the convolution of the thermal and instrumental profiles considered in each case, 
and the thick lines are the result of the convolution of the above mentioned lines, leading 
to the ’’synthetic observed profiles”. The striking (albeit obvious) fact illustrated in this 
figure is that intrinsic profiles narrower than Vf ^ 2 gg are gaussianized by the thermal and 
instrumental effects. Clear departures from a Gaussian profile are only seen for u/ > 3 
Gg- Note also that the values of Gear are systematically smaller than the intrinsic width, 
defined hj Vf. This is due to the convolution algorithm which by conserving flux leads, 
when considering a Gaussian and a steep profile (such as the law), to a re-distribution 
of the flux towards the base of the profiles. 

The figure allows for a direct assessment of the properties of the supersonic flow found 
in giant HII regions from the constraints that one can derive from the present study. Take 
for example the “quiescent flow” of 30 Doradus found as a Gaussian background component 
{gsodot ^ 18 kms“^) both in the core and the halo of the region (see CK), or the Geo 4 ^ 17 
kms^^arising from the main emitting knots in NGC 604 (see Sabalisck et al. 1995). Given 
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the expected thermal broadening in both regions 10 kms and the similar instrumental 
broadening (ainst 10 kms in both studies one can derive a (Tg ^ 14 kms Thus the 
supersonic flow in the background of both regions can be well explained with the results 
shown in flgure 6, in the panel with coordinates (ac = 10 kms^^, Vf = 20 kms^^). However, 
such a successful match immediately leads one to wonder why the terminal velocity (vf) of 
the unresolved shells is supersonic in such regions, a fact that implies a supersonic random 
speed of motions of their local interstellar medium. Note that ii Vf was smaller, say 10 
kms“^, the resultant profile would be narrower, and therefore different to those observed in 
our nearest giant HII regions. The situation becomes even more critical if one considers the 
full range of a values found in giant HII regions, which extends up to 40 kms^^. These, 
if produced in each case by a collection of unresolved wind-driven shells, takes us in our 
Figure 6 to even larger values of Vf and thus to faster supersonic random speed of motions 
of the local ISM in these sources to account for the missing slower shells. Furthermore, if 
the thermal and instrumental broadening remain ^ 10 kms^^, larger v/s also lead to wider 
flat-topped synthetic line profiles, in clear disagreement with the observations of giant HII 
regions which show a clear Gaussian behavior of the integrated emission lines across the 
whole range of supersonic a values (see Arsenault & Roy 1988). To gaussianize those line 
profiles it would be necessary to move towards the upper right-hand corner of Figure 6. 
This however, can only be achieved by assuming large changes in the thermal broadening, 
not only caused by a possible larger HII region temperature but rather due to the intrinsic 
supersonic random motions in each of the giant HII region sources. This brings us, once 
again, to infer the existence of another mechanism, different from the winds from massive 
stars, necessary to stir the ISM locally and thus explain the early disruption of unresolved 
wind-driven shells caused by massive stars while leading to a Gaussian integrated profile. 


4. Results and discussion 
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Given the large number of massive stars causing the photoionization of giant HII regions 
and their large mechanical energy output, through strong stellar winds and supernova 
explosions, the volume of a nebula mechanically unperturbed should become largely reduced 
as a function of time, and thus the observed line profiles should evolve to be dominated 
eventually by the emission from matter affected by the disruptive action from massive stars. 
Under these conditions the unperturbed gas may be rescued only by means of high spatial 
and spectral resolution, and in evolved HII regions high sensitivity instrumentation would 
be particularly required. All objects subjected to such studies (30 Doradus, NGC 604, NGG 
5471, NGG 5461, etc; see Munoz-Tunon 1994) show Gaussian profiles with a width, a a 
value, that implies a supersonic velocity dispersion in regions where the presence of strong 
stellar winds is not evident. We have explored the possibility that unresolved wind-driven 
shells produced by massive stars could account for the range of supersonic a values found 
in giant HII regions (15-40 kms^^). The main results from our analysis are : 

1) A collection of spatially unresolved, concurrently evolving wind-driven shells 
immersed in an HII region cannot cause a supersonic Gaussian line profile unless the 
wind-driven sources have a strongly peaked and thus peculiar age distribution (see Figure 
3). This can also be interpreted as a short-live event and thus, given the presence of 
supersonic line widths in all giant HII regions, makes this possibility rather unlikely. 

2 ) A uniform age distribution of the wind-driven sources leads to distinct non-Gaussian 
intrinsic line profiles and thus do not match the observations. These present a fiat-topped 
core defined by the range of terminal shell velocities {i.e. from -Vf to Vf), and a steep decay 
as up to the largest detectable shell speed (vi). 

3) Spatial undetectability of wind-driven shells is clearly favored by large background 
densities and low energy winds. For the closest giant HII regions, such as 30 Doradus where 
the supersonic gas is found both in the core of the region (with typical densities ^ 10^ cm^^) 
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and in the halo (where the typical density is ^ 1 cm^^), to have a collection of unresolved 
shells causing the supersonic width of the observed lines implies rather low-energy winds 
with properties similar to those typical of low-mass stars (E^ind ^ erg 

4) At the same time, the conditions that produce an intrinsic supersonic Gaussian 
profile while accounting for the unresolved shells as well as for a large range of values of wind 
energetics and background densities (e/n), lead unavoidably to large final shell velocities 
(vf). These imply supersonic values of the random speed of motions of the local interstellar 
medium disrupting the evolving shells long before they are capable of slowing down to 
speeds below 10 kms^b To our knowledge, the only possibility of causing and maintaining 
a supersonic value of the random speed of motions is through the continuous passage of a 
collection of low-mass wind-driven stars and their cometary shocks ( Tenorio-Tagle et al. 
1993). In such a case the motion of the ensemble is directly related to the total mass of the 
system, while the energy causing the stirring comes from the winds of the low-mass stars. 

5) Thermal and instrumental broadening tend to gaussianize line profiles, and the 
effect is more noticeable the narrower the intrinsic profiles are with respect to the Gaussian 
filter ((Jc). Note however, that the convolution of a narrow line arising from a system 

of unresolved shells and a Gaussian filter may lead to a broad Gaussian profile with a 
supersonic a value only if the final speed of shells (that defines the width of the fiat-topped 
core of the intrinsic profile) is also supersonic. This implies that the local random speed of 
motions of the target giant HII region is also supersonic, and thus capable of disrupting 
the evolving shells when Vg = Vf long before they reach, as in the Galaxy, a speed of < 10 
kms^^. Broader intrinsic lines (with vj > 30 kms^^) would however appear as described in 
2), and the only way to make such profiles acquire a Gaussian shape is through a convolution 
with a broader filter. The expected thermal and instrumental widths are, however, well 
determined and thus, the only remaining possibility is to add to the Gaussian filter a 
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supersonic component representing the local random speed of motions of the observed HII 
region. 

Thns onr attempt to explain the snpersonic Ganssian lines observed in giant HII 
regions with a large collection of unresolved wind-driven shells from massive stars leads 
in all cases to the existence of an underlying supersonic turbulent motion, independent of 
the wind-driven shells cansed by massive stars. Fnrthermore, this nnderlying snpersonic 
motion is what limits and defines the different terminal speed of shells in each of the giant 
HII regions. In our approach we had assumed static wind driven sources developing the 
nnresolved shells, althongh clearly an alternative way to ganssianize a broad emission with 
an intrinsic behavior is to allow for the motion of the sonrces within the gravitational 
potential of the whole (stars -|- gas) system. However, given the energetics of stellar winds 
and the densities of giant HII regions, only in the case of low-mass stars will the stellar 
wind confignrations be transformed into an ensemble of bow, or cometary, shocks raming 
through the nebula. These, while stirring the photoionized gas causing a supersonic random 
speed of motions, lead to the observed broad Ganssian line profiles. 

A fnrther indication that winds from low-mass stars are acting over large volnmes 
comes from the often disregarded fact that the HI clouds associated with giant HII regions 
also show a comparable snpersonic width, and both 30 Doradns and NGG 604 are good 
examples of this. McGee & Milton (1966) determined a a valne of 15.8 kms^^for their 
position (L32) which encompasses the massive (1.4 xlO^ Mq ), 500pc diameter HI cloud 
associated to 30 Doradus. There are also several other regions (positions L43, 48, 49, 50) 
with a snpersonic velocity dispersion all located in the neighboring kpc of 30 Doradns. 
Position L32 is however the most intense of these lines by more than a factor of 10. In the 
case of NGG 604, Wright (1971) measured a a of 23 kms“\ also largely supersonic for an HI 
clond. There are many other examples. Viallefond et al. (1981) determined and compared 
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the HI and HII supersonic a values of the giant regions in MlOl (NGC 5471, NGC 5462, 
NGG 5461), Hodge 58, NGG 5455, and NGC 5447, and all turn up to be rather similar. 

The measured ghi implies a large amount of kinetic energy in the HI phase, energy that 
cannot be attributed to the energetics from massive stars. The close agreement between 
the HI a values and those derived for the associated HII region gas, indicate therefore a 
different source of energy. 

All of these facts, i.e. the favoured low values of e, the large final shell velocities or 
the implied supersonic random speed of motions of the local ISM, and the similar values 
of a found for the ionized and neutral components of each giant HII region, are elements 
contemplated in the cometary stirring model proposed by Tenorio-Tagle , Munoz-Tunon & 
Cox (1993). The model predicts a gas ^ cr stars, as has been confirmed for the Orion cluster 
and the Orion molecular cloud (see review of Genzel & Stutzki 1989). For giant HII regions 
however, this fact still awaits ratification. 
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Fig. 2.— Unresolved wind-driven shells. The emission line profile resulting from a collection 
of photoionized unresolved wind-driven shells with a uniform age distribution and initial and 
final speeds of 100 and 10 kms^^, respectively. Note that larger values of the initial speed 
would hardly affect the resultant profile. 
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Fig. 4.— Wind-driven shells. Total intrinsic line profiles resulting for a system of shells with 
different terminal velocities between (10 kms“^) and (30 and 50 kms“^), and a 
power-law weight function. 





Fig. 5.— Wind-driven shells. Minimnm shell velocity (vs) plotted against the ratio of e/n (in 
units of 10^® erg s“^). Shells become detectable as they exceed the Vcrit value determined by 
the solid lines and cross into the area below the curves. The range of e/n (assuming an Vcrit 
= 10 pc) for the detected slow expanding shells detected in 30 Doradus (30 < 50 kms^^) is 
indicated by the dashed area, while unresolved shells covering the range of a values detected 
in giant HII regions (15-40 kms“^) is indicated by the dotted area. Note that the latter 
becomes smaller, causing extreme values of e/n as the spatial resolution increases {i.e. for 
smaller values of Vcrit) 
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Fig. 6.— Synthetic line profiles. Predicted v~^ line profiles (narrow lines) for a multiple- 
shell system convolved with a thermal plus instrumental Gaussian (dotted lines) of different 
widths. The thick lines show the resulting normalized line profile after convolution. The 
corrected values of a are indicated in each frame. 























